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Abstract 


The discharge characteristics of a Li/Li,CoO> cell were simulated by numerical calculations. Based on the program proposed by Newman’s 
group, the change in entropy due to lithium-ion insertion into the active materials was introduced into the open-circuit potential (OCP) for more 
practical applications. The OCP of LiCoO, (y= 1) at a full-discharge state was 3.82 V, which was determined from the original program. However, 
it should show a rapid decrease to minus infinity when discharge is complete (y= 1), which is due to the introduced entropy term. In addition, the 
local use of active materials across the Li CoO, electrode was not uniform, whereas the original program showed a constant value throughout 
the electrode. Based on the present results, the introduction of an entropy term to the OCP of Li,CoO, works reasonably well for more practical 


applications. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Rechargeable lithium-ion batteries have been used in a wide 
variety of portable electric devices due to their high energy den- 
sities. However, these batteries still need to be improved before 
they can be used in hybrid electric vehicles (HEV) [1,2]. Rapid 
charge and discharge reactions are required when lithium-ion 
batteries are used for high power applications. Thus, lithium-ion 
transfer must be very rapid in such batteries. Lithium-ion transfer 
in lithium batteries consists of four different steps; lithium-ion 
diffusion through positive electrode materials, lithium-ion trans- 
port through an electrolyte, and lithium-ion transfer across the 
interfaces between the electrolyte and active electrode materi- 
als. Of these steps, the diffusion of lithium-ion through active 
electrode materials is well known to be slow. In addition, the 
charge transfer resistance at the electrode/electrolyte interface 
is expected to be low by extending the effective surface area of 
the active materials. Therefore, fine particles of battery active 
materials (ca. 5—10 um) are used to reduce the diffusion path in 
commercial lithium-ion batteries [3-5]. In lithium-ion batteries 
intended for high power use, the particle sizes of active materi- 
als still need to be reduced further [6]. We previously prepared 
nano-sized LiCoOy particles by an excess lithium method or 


* Corresponding author. Tel.: +81 92 583 7790; fax: +81 92 583 7790. 
E-mail address: yamaki@cm.kyushu-u.ac.jp (J.-I. Yamaki). 


0378-7753/$ — see front matter © 2007 Elsevier B.V. All rights reserved. 
doi: 10.1016/j.jpowsour.2007.06.089 


by laser ablation, and reported that the rate performance was 
improved by use of the resultant nano-sized LiCoO2 compared 
with that using micron-sized LiCoO2 powder [7,8]. However, 
the cell performance was influenced not only by the particle size 
of the active materials, but also by other geometric factors such 
as the thickness of the electrode. Such inherent complications 
should be mainly due to the porous structure of the electrodes; 
the ohmic potential drop and mass transfer were very intricate, 
and there was no way to separate them. Some powerful meth- 
ods are needed to review the experimental results methodically, 
which should contribute to the design of high-power batteries. It 
may be useful to simulate cell performance by calculation. Doyle 
et al. proposed a mathematical model of a Li/insertion electrode 
cell using a porous electrode theory, and reported the simulation 
results [9,10]. Newman et al. studied the charge/discharge char- 
acteristics of a dual-insertion cell of C/LiMn2O4 by numerical 
calculation, and compared them with experimental data [11,12]. 

In the present study, the discharge characteristics of a 
Li/Liy,CoO 2 cell were studied by numerical calculation using 
the program proposed by Newman et al. [13]. Based on that pro- 
gram, we introduced changes in entropy to the OCP for more 
practical applications. 


2. Model development 


The program used (dual. f (version 4.0)) is freely available 
from Newman’s group through their website. The program is 
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based on a cell model consisting of a porous positive-electrode, 
a separator, and a lithium metal negative-electrode. The porous 
electrode was composed of sphere-shaped particles of active 
materials, inert conductive additives, and polymer binder. The 
void space in the electrode was filled with an electrolyte solu- 
tion. In the model, the one-dimensional transport of lithium-ion 
from the negative electrode to the positive electrode through the 
separator was considered, and the electrochemical insertion of 
lithium-ion was assumed to occur at all points on the active mate- 
rials. The program was composed of six differential equations, 
which are described elsewhere in detail [11,13]: 


(1) material balance on the electrolyte; 

(2) current flow through the electrolyte; 

(3) current flow through the active materials; 

(4) pore-wall flux subjected to the Butler—Volmer equation, 
which is related to the divergence of the current flow in 
the electrolyte; 

(5) Fick’s first law for spherical diffusion; 

(6) Fick’s second law for spherical diffusion. 


These equations were solved simultaneously to obtain the 
following solutions: 


(1) lithium-ion concentration in the electrolyte, 

(2) lithium-ion potential in the electrolyte, 

(3) lithium-ion concentration on the active material particles, 

(4) ionic current flowing through the electrolyte, 

(5) pore-wall flux across the interface between the active mate- 
rials and electrolyte, 

(6) surface potential of the active materials. 


In the present work, LiCoO2 was used as the active mate- 
rial for the porous positive-electrode. Accordingly, the overall 
reaction can be described as follows: 


LiyCoO2 + (1 — y)Li = LiCoO2 


The cell model is shown schematically in Fig. 1. The thickness of 
the negative electrode, the separator, and the positive electrode 


Cu current collector (25 um 


Limetal negative-electrode 


o 
separator (25um) 
25 
LiCoO, porous positive- 
electrode (100 um) 
125 
Al current collector (25 um) 
x (um) 


Fig. 1. A schematic model of a lithium battery. 


were set to 100, 25, and 100 um, respectively. The surface of 
the negative electrode was defined as x= 0. The diffusion coeffi- 
cient of lithium-ion within LiCoO) particles and a standard rate 
constant (k°) were determined based on the literature [14]. The 
former was assumed to be constant at 1.0 x 10~!© m?s~!. The 
latter was evaluated to be 1.5 x 107° A mol?” m°” from the 
value for charge transfer resistance (Rot) (Eq. (1)): 


R RT (1) 
ct Fio 

where R is a gas constant, T the absolute temperature, and F is 

the Faraday constant. 

In this study, the effects of particle size of LiCoO2, porosity 
in the porous positive-electrode, and conductivity of the positive 
electrode on the cell performance were investigated. Hence, the 
preset parameters in the original program were used except for 
those mentioned above. The change in the electrode potential 
of the LiyCoOz electrode with lithium content (y), which was 
obtained from the original program, is shown with a broken line 
in Fig. 2. The OCP decreased from 4.7 V toward 3.8 V with an 
increase in y, but never fell below 3.82 V. In addition, the elec- 
trode potential varied linearly at around y=1 with an increase 
in y, whereas a rapid decrease in electrode potential is seen in a 
real cell [15]. The Butler-Volmer equation (Eq. (2)) used in the 
program was 


oo, 0.5F(n — U) —0.5F(n — U) 
j= jo {ex | exp |} (2) 


RT RT 


where j is the current density, 7 is an overpotential defined as 
n=, — D2 (Pj, D2: electrical potential on the active material 
surface and in the electrolyte, respectively), and U is the OCP. 
An exchange current density of jo is defined as 
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Fig. 2. Variation of the open-circuit potential of the LiyCoO2 electrode before 
and after changes in entropy were introduced into the OCP of Li,CoOd. Particle 
size of LiyCoO2: 100nm, conductivity of the positive electrode: 100S m7!, 
porosity in the positive electrode: 0.772. 


T. Doi et al. / Journal of Power Sources 174 (2007) 779-783 781 


where Csmax is the maximum concentration of lithium-ion in 
LiCoO? (y= 1), Css is the lithium-ion concentration on the sur- 
face of LiCoO} particles, and C is the lithium-ion concentration 
in the electrolyte. In the present work, the change in entropy 
due to lithium-ion insertion into the active materials was intro- 
duced to U as follows [16,17]; as a first order approximation, 
interactions among Li-ions in LiyCoO2 were ignored: 


1 RT H 
U' = U — — In | —— (4) 
F 1-y 


U’ was then substituted for U in Eq. (2). 


3. Results and discussion 


The change in the OCP with y after introduction of the entropy 
term is shown with a solid line in Fig. 2. The electrode potential 
showed a rapid decrease at around y= 1 and fell below 3.82 V, 
which is very similar to the results obtained with a real cell [15]. 
The entropy term also varied rapidly at around y = 1, as shown in 
Fig. 2. These results indicate that the OCP curves were affected 
by the variation of entropy with y. 

Fig. 3a shows the variation in discharge capacity (broken line) 
with current density, as obtained with the original program. Cell 
voltage, positive-electrode and negative-electrode potentials at 
the end of discharge are also shown with solid line in Fig. 3a. 
The cell voltage was defined as the difference between the elec- 
trode potentials. The particle size of LiCoO2, conductivity of the 
positive electrode, and porosity in the positive electrode were 
set to 100nm, 100 S m~}, and 0.772, respectively. The positive- 
electrode potential at the end of discharge fell to only 3.82 V 
even at a low current density of 1.0mA em"). and remained 
above 3.5 V at current densities of less than 15 mA om". On the 
other hand, the negative-electrode potentials increased with an 
increase in the current density, which should be due to overpoten- 
tials. As a result, the cell voltage never reached 3.5 V at current 
densities of less than 15 mA cm". Fig. 3b shows the variation in 
discharge capacity, cell voltage, positive-electrode and negative- 
electrode potentials with current density that was obtained with 
the program after the entropy term was introduced into U. Even 
at a low current density of 1.0 mA cm~?, the positive-electrode 
potential dropped to 3.49 V at the end of discharge and the cell 
voltage fell to 3.5 V. Although the potentials of both electrodes 
increased with an increase in the current density, the cell volt- 
age still dropped to 3.5 V at the end of discharge regardless of 
the value of the current density; since discharge was cut off at a 
potential difference of 3.5 V and the negative-electrode potential 
increased with increasing current density, the positive-electrode 
potentials at the end of discharge increased with an increase in 
the current density, as shown in Fig. 3b. A discharge capacity 
of 137 mAh g7! was obtained at 1.0mAcm~?, and gradually 
decreased to 130 mAh g7! at 15 mA cm~? followed by a signif- 
icant decrease at current densities above 20 mA em". 

Fig. 4 shows the use of active material in the Li,CoQ? elec- 
trode at the end of discharge. The original program gave a 
constant value of y = 1.0 throughout the entire electrode (broken 
line). On the other hand, the introduction of the entropy term 
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Fig. 3. Variation of discharge capacity and the cell voltage at the end of discharge 
with current density (a) before and (b) after variation of entropies was introduced 
into the OCP of LiyCoOp. Particle size of LiyCoO2: 100 nm, conductivity of 
the positive electrode: 100 S m~!, volume fraction of electrolyte in the positive 
electrode: 0.772. 
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Fig. 4. Local use of active material in the LiyCoOz electrode at the end of the 
discharge reaction. Particle size of LijCoO2: 100 nm, conductivity of the positive 
electrode: 100 S m~!, porosity in the positive electrode: 0.772. 
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Fig. 5. Changes in discharge capacity with current density at various porosities 
in the positive electrode (ranging from 0.318 to 0.772). The particle size of 
LiyCoO2 was (a) 10 nm, (b) 100 nm, and (c) 1 um. Conductivity of the positive 
electrode was 1.00 S m~!. 


resulted in a difference in y across the electrode (solid line); y 
gave a smaller value of 0.87 in the middle of the electrode than at 
both ends of the electrode, which should be due to the influence 
of ionic and electric conductivity. These results provide evidence 
that the introduced entropy should work reasonably well in the 
program to affect cell performance. These results indicate that 
the introduction of an entropy term into U’ makes it possible to 
control the discharge reaction within a range of predetermined 
voltages. 

The influence of porosity in the positive electrode on cell per- 
formance was investigated by using the program that included 
an entropy term. Fig. 5 shows the variation in discharge capac- 
ity with current density at various porosity values ranging from 
0.318 to 0.772. The particle sizes of LiCoO? used in Fig. 5(a)-(c) 
were 10nm, 100nm, and 1 um, respectively. The discharge 
capacity tended to increase with a decrease in particle size. The 
porosity at which a maximum capacity could be obtained dif- 
fered depending on the particle size of LiCoOz; the particle size 
of 10 nm gave the largest capacities when the porosity was 0.545, 
and the corresponding values were 0.431 and 0.318 for 100 nm 
and 1 um, respectively. These results suggest that there should 
be an optimal value for the porosity of a porous electrode to 
enhance the rate performance of batteries. 


4. Conclusion 


Based on the program proposed by Newman’s group, the dis- 
charge characteristics of a Li/LiyCoO cell were simulated by 
numerical calculations. Changes in entropy during the insertion 
of lithium-ion were introduced into the OCP, and their effects on 
the discharge characteristics were studied. The discharge volt- 
age fell to 3.5 V, and a rapid decrease in the electrode potential 
was seen at around y=1 with an increase in y by introducing 
the entropy term, which are similar to the results obtained in a 
real cell. In addition, nonuniformity in the local use of active 
materials across the LiyCoO2 electrode was observed. These 
results indicate that the entropy term we introduced should work 
reasonably well for more practical applications. The discharge 
performance was investigated using LiCoOy particles of various 
sizes, ranging from 10 nm to 1 wm. The porosity that gave the 
largest capacity was different depending on the particle size of 
the active materials. 
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